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12,  Abstract;  The  fundamental  equations  of  nngnetofluiddynamlcs  are  derived 
For  incompressible  fluid ,  itia gnetobydrodynandcs ,  the  Important  parameters  are  the  Heynolds 
number  the  magnetic  pressure  number  which  is  the  ratio  of  magnetic  preeaure  to 
dvmmic  preasu’-e,  and  the  velocity  number  R^,  which  1?  the  ratio  of  the  fluid  velocity  to 
the  cK-'rocterlstic  velocity  of  which  the  magnetic  field  is  moving  through  a  conductor. 

Sciss  exact  solutions  and  properties  of  the  equations  of  fnagnotohydrodynamlcs  are  given. 
Stability  of  ?.r.mlnar  flow  and  turbulence  In  raafjrie to f^pr> dynamics  are  briefljr  reviewed^. 
Finally  aotrrt  rwignoto hydrodynamic  experiments  are  deacribedo 

For  ccnprcsnlblo  fluid,  magneto gasdynaraics,  the  impoitant  parameters  are  still  R^, 
and  Ry  plus  other  well  known  In^wrtanC  paranetero  of  ordinary  gas  dynamics  such  as 
Mach  number,  Prandtl  number  and  ratio  of  specific  heatso  Both  the  waves  of  smell 
amp'iitjlo.  Alfvtm's  waves  in  compressible  fluid,  and  shock  unveo  in  nvtrrictogasdynamic* 
arc  di;tcui3E;';d ,  , 


Best  Available  Copy 


Intr^sductloao 

Electromaf  latic  phiMif  rtanfi  in  solid  coriduc:tDis  have  bean 
well  known  for  a  lonp:  1101(0  but  th3  e].ectrct:naf natic  phanoraeoa 
in  fluids,  llcuid  ov  gases,  are  not  so  well  known^  Only 
recently  the  problain  of  rnagnetohydrodynanics  has  been  attract- 
ing  the  attention  of  some  reseai-ch  workers  because  it  is 
irpo riant  in  astrophysics,  geophysics  and  the  behavior  of 
interstellar  gas  masiesc^*^  The  main  diffeitnco  in  the  electro- 
m''gnetlc:  phen.-.:'.  na  io  flav.ls  frcm  those  ir.  sc  l  lcis  i;.,.  due  to  the 
fact  that  necha  lical  forces  dei'lv'ng  f  rt/in  elf-otrlo  currents  nay 
produce  f?.uid  d  riamlo  motions  arJ  the  x’luld  dy carle  motions  nay 
produce  electro  lagc&tic  pheronor  >,  In  other  wcixls.,  tiiere  Is  on 
interaction  bet  .een  the  electroir  ■  nstto  forces  rr\d  the  ordinary 
fluid  dynamic  f.-roco  Tnis  inter  ctlon  phenc'mona  Vicoraee 
irportant  whonerer  tJie  electromsg-atlc  forces  are  of  the  sane 
ordar  of  magritode  as  the  inei'ti;  1  forces  or  viscous  forces  of 
the  fluids  Actoally  the  discovery  of  nsagnetohy'drodynaniica  was 
made  in  cosmic  physicss.  Magnetic  phenomena  of  different  types 
has  been  observed  and  needed  an  explanations  In  ordinary  labor» 
atory  experlraente  of  dis^targss  of  electrical  current  in  gases, 
the  nschanical  effects  are  usually  small,  the  results  of  these 
experiments  cannot  be  used  to  explain  the  electromagnetic  phen¬ 
omena  affecting  the  notion  of  gaseous  masses  of  cosmical  dimen¬ 
sions.  Hence  study  of  electromagnetic  phenomena  in  liquid  con¬ 
ductors  has  been  extensively  carried  out  where  In  some  respects 


the  conditions  are  more  sli.ulrr  to  those  occuri'ing  in  cosmical 
physic® »  Thuj  che  magriotohydi'odynamlcs  has  been  fo’antU  The 
result  of  mafpietohyrJi-odynanics  has  also  gottb  ongireerlnK 
interest  because  of  tie  utility  o:C  induction  flow-meters,  which 
rely  on  the  generation  of  a  measurable  potential  difference  in 
the  fluid  in  a  iirection  perpendicular  to  the  motion  and  to  the 
magnetic  field* 

Recently  tacausa  the  interest  of  hypersonic  flow  of  missiles, 
the  interaction  of  the  alectromagnatic  force  and  fluid  dynamic 
forces  in  gases  are  not  negligible  and  may  be  pit>duced  in  lab- 
oiMtories  by  the  use  of  shock  tube.  Ihus  aom  preliminary 
studies  of  ma  ^netog^jsdynaaics  has  been  initiated*  9ut  nuch 
has  to  be  done  -cforo  a  conplete  i  nderstanding  of  these  phen¬ 
omena  is  possible.  2h  this  pressr  t  paj^er,  I  shall  make  a  brief 
reTlev  of  the  pnesent  status  of  mfignetohydrodynanios  and  nagneto- 
gas dynamics a 

2,  P^indamental  Equaxions. 

In  studying  the  magnetic  fluid  dynamics,  one  has  to  deal 
with  the  fluid  (iynamical  equations  and  the  electrodynamic  eque- 
tlons  sifflultaneuuely.  In  the  general  three  dimensional  flow 
there  are  16  unlmcvne  in  the  problems  of  a  magnetofluiddynanlcs, 
i.e.  j 

(a)  Ihe  magnetic  field  strength  ^  (3  components) 

(b)  The  electric  field  strength  E  (3  ooiiponents} 

(c)  The  electric  current  density  J  (3  components) 


(d)  'Hie  cxc33n  .jloctrio  charga  {% 

(e)  Hie  fliid  velocity  vector^  (3  conponents) 

(f )  The  pressure  of  the  fluid  ^ 

(g)  Ihe  density  of  the  fluid  p 

(h)  The  temperature  of  the  fluid  To 

Fbr  these  unknowns,  we  ha  we  to  find  16  relations  which  are 
the  fundamental  equations  of  magnetoflulddynamlcc  and  which  are 
glvan  below t 

The  Ifexweir  equations^  ars  atsuned  to  be  true  in  raagneto- 
fluiddynanicso  'Itoj'  .'iro  six  In  n'lt'.ber  and  are  as  follows  t 

(2.1) 


(2.2) 
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J 


£  - 


where  ^  is  the  gradient  operatoro  Gibbs  vector  notatione  are 

usedo  t  is  the  timeo  €  is  the  dielectric  constant  and^^ia 

t\KS 

the  magnetic  permeabllityo  The  jfSR.  uxiit  system  Is  usedo 
The  current  density  equatipn  is 


r,.'  •'  ■ 


(  (2.3) 

where  Q  is  the  electric  conductivity,  The  terms  with  q  aze  the 
coupling  terms  ulth  the  fluid  dynamic  equations  which  represent 
the  interaction  phenomena  <> 


The  conservation  of  electidLc  charge  gives 


^  (7 
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(2.U) 
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Equations  '  2.-1)  to  (2„li)  are  the  electrona.^etic  equations 


vi  th  coupling  torrrts  .iuc  to  fluid  djTiamic  motiono 

The  conEer<. ation  of  mans*  i,e«,  the  equation  of  continuity 


jji  fluid  dynamics,  la 


.  -  .  —  ^  .  .  -r  V  s  ( 

-/  •  * 

Fbr  incompressible  fluid, f  :  constant  ^  /q  t  equation  (2,5) 
becomes 


(2<.5) 


-  J 


(2«5«) 


The  equation  of  motion  is  as  roUowa 


where  *is  the  stress  tonaor  due  io  vincosity  which  is  usually 
assumed  to  be  linear  homogeneous  function  of  rate  of  change  of 
velocity.  The  1  3  co^iponent  of  th:  stress  tensor  C  *i»y  b« 
written  as  foUous: 


t;:  ^ f^f  ^  s’Ja'  '.-  j* 


(2»7) 


where  1,J  s  l,2,3o  Subscript  i  or  ^  refers  to  the  1  or  j  oojik 
ponent  of  the  vector  respeotlvelyo  Is  the  ith  oonponent  of 
velocity  vector  q  and  is  the  ith  component  of  cartesian  spatial 
eoordinatso  s  0  If  i  j,  1  if  i  »  jo  /'  is  the  ordinasy 

coefficient  of  vlscosityo  We  assure  that  the  second  coefficient 
of  viscosity  is  For  incoiqpressible  fluid  with  constant 

viscotilty,  the  viscosity  tame  uey  be  greatly  8iiq>lifiad  mad 
equation  (2.6)  becomes 


(2c  6«) 


t :  -  (f  Jx  ^  “f 

whemV-^s^’V. 

It  is  an  orapirlcal  fact  that  thore  is  a  runctional  relation 
between  the  density  f  ,  the  pressure ^  and  tiio  tejape nature  T 
of  a  fluids  This  relation  is  known  as  equation  of  states  fbr 
gas  dynamics,  the  perfect  gas  lav  is  usually  used  wlxich  is 

J^fRT  (2.e) 

whore  R  is  the  gas  constant.  In  magneto  gas  d3na!tilcs  equation 
(2,8)  is  also  used  as  the  equation  of  state. 

In  nagnetohydrodynamics,  eqiiation  (2,3)  is  replaced  by 

)  "  “  constant  (2,8a) 

lAlch  is  assumed  to  be  known. 

The  last  relation  is  the  energy  equation  vhich  is 


-f*  -  i£f  —  g  "  j .  o 


(2.9) 


where  the  sonmatlon  convention  is  used,  i.e. 


^  z  ^ 


and 


5  |-If-0E-^% 

I  =  Internal  energy  of  the  fluid  per  unit  mass 
s  C^T  for  perfect  gas,  *  spc-cific  heat  at  constant  volune, 
I]^  5  ■—  ^l®i  “  electric  energy  per  unit  volume. 


D,,  =  I  mai^istic  eiforgj  per  unit  wla-^o 

c^<? 

irate  of  energy  produced  by  external  agencieu.  For  adiabatic 
system  s 

z  Ij  conponant  of  viscous  stress  tensor  givan  by  (2o7)o 

2  ij  coittponant  of  electrora gna tic  stress  tenst^r 

.  e  £y  V  H,-  -  -tie  ' 

K  S  coefficient  of  heat  conductivity 

s  ith  ooniponant  of  the  Rjyrtlng  vector  Sj  where 
?  X  T 

component  of  the  x'adlal^on  energy  flux© 

%  -  Ra 

s  rate  of  radiation  energy  emission  per  v.nit  volume* 

%  S  rate  of  radiation  energy  absorption  per  unit  volume* 

3h  the  case  of  Incompressible  fluid,  l*e*  magnetohydrodyna- 
rales «  if  ve  assume  that  the  coefficient  of  viscosity y*<  is 
constant,  we  may  solve  the  unknowns  Hj  /  /e  ,  «nd^ 

from  equation  (2,1)  to  (2*6).  without  considering  the  energy 
equation.  After  these  'unknowns  are  found,  equation  (2.9)  gives 
the  temperature*  For  conpressible  fluids,  i.e*  magneto^* 
dynamics,  ve  have  to  deal  with  equations  (2,1)  to  (2,9)  slimil- 
tansously*  Li  the  following  we  shall  review  briefly  the  present 
status  of  nagnetohydrodynamics  and  magneto gasdynaalcso 


Fart  lo  MAGUETO-^HYDaiDv-MiriCS 


Bnportant  Paramaters  In  r6i{?ietohydrodyiirmtC8o 
Vlu  cons  id -^r  only  the  cases  nhere  the  velocity  of  the  flow 
"q  is  much  s/nal  jer  tJian  the  velocity  of  li.Tht  Co  The  relati¬ 
vistic  effects  Jnay  be  omitted.  Since  the  cr.ergy  in  the  electric 
field  is  of  th:;  order  of  of  energy  in  the  magnetic 

field,  iii  this  ease,  the  energy  in  the  electric  field  can  be 
neglectedo  Cc  seqnently  in  ordinary  nagnetohydrodynaraics,  we 
considar  only  he  interaction  between  the  velocity  field  q,  and 
the  magnetic  f  eld  h,  Ih  this  approximation  we  may  put  A  s  0 
In  the  flmdamantal  equations* 

The  fundaiiental  equations  for  magneto  hydrodynamics  nay  be 
obtained  fzom  equations  (2.1)  to  (2,6)  in  the  following  non- 


di'cenalonal  f  o 


(3*1) 


(3.2) 


(3.3) 


luiiere  the  non-dimensional  quantities  are  defined  as  follows  i 


r-fu-  > 


--tL 


U  is  the  characteristic  velocity,  L  is  the  characteristic 


length  and  is  tho  charactorrlstlc  magnetic  field  of  this 
system.  Hie  probleiri  is  to  solve  tho  unknowns  p*,  T*  ani'^ 
in  terras  of  the  important  paraTneters  R^. 

The  meaningB  of  these  throe  Important  parameters  are  as 
follows : 

(1)  is  the  well  known  ^teynolds  nuniser  which  character¬ 


ises  the  viscous  flow. 

f  UL 

n  s  Inertial  force  s  . —  (3.5) 

Wbeo'5"7o~?5r  ^ 

(2)  is  the  magnetic  pressure  number  which  is  the  ratio 
of  magnetic  pressure  ~  *'a  dynamic 

pressure  11^/2. 


-  ygnetic  pressure  -  ^ 


lynamle  pressure 


4  (/"A 


(3.6} 


Only  idien  is  of  the  order  of  unity  or  larger  the  fluid  flow 
will  be  affected  by  the  magnetic  field.  If  ^  <  1^  tbs  terns 
due  to  magnetic  field  in  the  equations  of  motion  may  be  ns^etsd 
and  the  ^vUl  not  be  affected  noticeably  by  ^ 

(3)  is  the  velocity  number  which  is 


/?,=  «  S- 


(3.7) 


where  =  f/TH  *  ^  regarded  as  a  charaeterlstle 

velocity  of  which  the  magnetic  field  is  moving  through  the  con¬ 
ductor.  If  U  >  >  the  field  is  practically  conqpelled  to 
follow  the  motion,  the  iragnetic  field  will  be  greatly  influsnoed 


by  th*  motion*  On  the  other  handl«  if  U  V^,  the  magietle 
field  will  not  be  influenced  noticeably  by  the  motion*  This 
is  the  zeason  idiy  the  raa gnetohydrodynamic  dl>?i8n3lon  L  is  veiy 
large  so  that  V^,  is  very  small* 

it*  Solutions  and  Properties  of  the  Equations  of  Hagnetohjydro- 
dynamics* 

It  Is  well  known  that  there  is  no  general  mthod  of  solu* 
tion  for  hydrodynamlcal  equations  alone*  Thus  it  is  much  more 
difficult  to  ln^^-estlgate  mgnetohydrodynamlcs*  One  of  the  first 
steps  to  understand  reagnetohydrodynamic  equations  is  to  find 
some  simple  exact  solutions  of  nagnetohydrodynamlcal  equations 
in  order  to  faring  out  its  essential  features*  One  way  to  find 
such  exact  solutions  will  be  the  generalisation  of  some  well 
known  results  of  ordinary  hydrodynamics*  The  following  gives 
some  of  these  results t 

(a)  TWo  dimensional  steady  flow  between  parallel  i^tee* 

Oonsider  two  parallel  plates  situated  in  the  plane 
X2*  s  d:  1  respectively,  being  one  of  the  non-dimensional 
spatial  coordinates  x^**  let  only  one  component  (x^*  -  vise) 
of  velocity  be  different  from  sero  and  be  a  function  of 
only,  i*a*  Uj^*(x2*)<»  FUrthernwre  we  assume  that  there  is  a 
uniform  extemaTi  magnetic  field  which  is  in  the  Xg*  -  direc¬ 
tion*  It  should  be  notice!  that  if  Is  in  the  Xj^*  or  x^*  - 
direction,  the  velocity  field  will  not  be  influenced  by  the 
magnetic  field  in  this  case*  Wilder  these  conditions,  the 


equations  (3ol)  to  (3.3)  glva 


c/  iy/- 
d  y.i 


:r  O 


(^•1) 


.  I  /V  *  I  *• 

f  ^  f?  P  R  ^  ^  <g~ 

where  ^  Ae  -»  - 

Ihe  general  solution  of  equation  (U.l)  la 

s. ^  CirJLAx^  ^  +  C 


(I1.2) 


(1*.3) 


i^re  A,  B  anl  C  are  arbitrary  constants  to  be  detenilned  fztm 
the  boundary  conditlon8« 

<1)  Plane  Oonette  flow.  The  boundary  conditions  are  X2*  m  0 
Uj*  •  Oj  *2*  *  Uj^*  a  1,  ^2*  s  -1,  uj^*  •  -I.  Equation 
(U.3}  becomes 

ui*  s  5^-|iE2*  (UJ») 

Ihla  flow  corresponds  to  the  flow  between  two  parallel  Infinite 
plates  novlnf;  in  oi^xisite  direction  with  velocity  i  0  and  at  the 
location  3^*  e  ±1  respectively.  There  Is  no  pressure  gradient 
in  -  direction. 

£f  "k”  tends  to  be  aero,  equation  {h»h)  becceies 

•  X2*  (^*5) 

whidh  Is  the  linear  velocity  distribution  of  ordinary  hydrod]^ 
nanios. 


If  "Id*  is  very  large,  ,  than 


-11. 


^  ^  -**=44^  — * « 

except  X2*-~>  i  lo  For  very  large  value  of  k,  in  the  central 
portion  between  the  two  plates »  the  velocity  Is  alnost  sero  and 
then  Increases  ^'ery  rapidly  to  the  values  ^  1  near  X2*  etl* 

Ihis  is  a  typical  phenomenon  In  magnetohydrodynaalcs  of  large  k« 
Tne  reason  Is  due  to  the  fact  that  the  magnetic  field  increases 
the  total  shearing  stress  of  the  flow  field.  Near  the  vall^  the 
shearing  stress  due  to  the  magnetic  field  tends  to  sero»  the 
total  shearing  stress  will  be  largely  produced  by  the  viscous 
force  only,  a  result  there  will  be  a  large  velocity  gradient 
near  the  wall.  This  effect  la  qualitatively  similar  to  that  of 
turbulent  flow.  In  turbulent  flow  the  total  shearing  etaess  of 
the  flow  is  much  larger  than  that  of  the  corresponding  laminar 
floe.  But  the  turbulent  stress  will  be  sero  at  the  vail;  tboe 
there  vlU  be  a  large  velocity  gradient  near  the  wall  of  turv 
bolent  flow  to  respond  to  the  large  shearing  stress  of  the  flow 
field. 

(ii)  Flans  PoiseuUle  Flow,  The  boundary  oonditlone  are 
*2*  «  il,  Uj*  «  Oj  Xg*  •  »1*  *  1*  Bqnatlon  (U.3)  boe<»as 


If  k-40,  (li.7)  becomes 
Uj*  s  1  •  *2*^ 


(lu7) 


(U.8) 


This  is  the  plans  PolseulUe  flow  of  ordinary  hydrodynamics.  Xf 
k-^eo. 


o 


exe«pt  X 


-4— 


±  1,  hance 

^  *  S'  1 


■» 


ih.9) 


exempt  x^  s  Xe  Here  again  va  have  a  large  velocity  gradlant 
near  the  valleo 

After  the  value  of  le  obtained^  the  induced  nagnetie 
field  in  x^*  *»  direction  and  the  preseure^  *  may  be  obtained 
by  quadratuz«« 


(Ik  ID) 


vhere  the  boundary  conditions  X2* 

«n  -  rWxf 


1 1,  :  0  are  used  and 


ehere  and  are  constants*  Is  the  pressure  gradient  in 
Xj^*  •  dlieotion,  (dp*/dai^*)*  ^  is  the  pressure  at  x^*  a  0, 

X2*  =:  £  lo  It  is  interesting  to  notice  that  x^*  wise  pressure 
variation  Is  the  sans  here  as  that  in  the  ordinary  hydrodynanics* 
The  Magnetic  field  inbroducw  the  variation  of  pressure  in  tbs 
X2*  •  direction*  this  point  again  is  slnilar  to  the  case  of 
turbulent  flow  of  ordinary  hydrodynamics* 

(b)  Other  simple  steady  flows*  There  are  many  other  slaple 


steady  flow  exact  solutions  had  been  found,  e.g*t 


(I)  Plov  5j)  a  circular  pipe  under  an  external  radial  nag- 
netle  field,  (ref,  8) 

(II)  Plow  In  a  rectangular  channel  under  transverae  nagaetle 
fields  (ref.  9)c 

(ill)  Plot!  between  concentric  circular  cylinders  under  oc 
stanb  axial  magnetic  field  (ref*  10)*  In  this  ease*  the  teloel.^ 
distribution  Is  not  affected  by  the  presence  of  the  magnetie 
field* 

(e)  Analogy  between  magratic  field  and  Tortlclty. 

There  is  a  formal  analogy  between  the  rortlelty  ^  of 

ordiiuQry  bydrodynanlcs  and  the  magietlc  field  H  of  nagnatohydrs- 

dynamics'^  0  The  fbndamental  equations  for  vortielty  O  art 

V'  a?  O 


^  -  V>i(fxZT)  =  ^ 


(lul2) 


The  fundamental  equations  for  the  magnetic  field  R 

y7>  ^  =  O 


dH 


3^  -vx(fAjr)= 


(U.13) 


irtiere  magnetic  dlffbsivity* 

Since  equations  (2ul2)  and  (lull)  ere  identical  in  form* 
we  may  i^ply  all  of  the  known  theorens  of  vortielty  in  ordinary 
hydrodynanlcs  to  the  magnetic  field  of  nagnetohydrodynamlcs* 


Fbr  instance*  from  Rslmholta's  theorem*  we  nay  show  that  vhen 
eonduetlrity  C  is  infinite*  i.e*^*^ :  0*  the  lines  of  magnetic 


force  move  irlth  the  fluid 


5*  stability  of  Laminar  Flow  In  MBgnetohydrodynMica« 

It  has  now  been  generally  recognised  tiiat  turbulent  notion 
la  the  more  natural  state  of  fluid  flov^  and  the  laminar  motion 
occurs  only  the  Reynolds  number  is  so  low  that  the  deyiatlone 
from  it  are  llnbls  to  be  damped  out*  In  cosmic  conditions* 
idiere  the  problem  of  magnetohydrodynamlcs  is  liq>ortant*  because 
of  large  dimensions*  l.e*  large  Reynolds  numbers*  one  vould 
expect  that  the  flow  will  be  mostly  turbulent*  Ihere  are  two 
problems  of  great  Interest  in  these  conditions*  One  Is  the 
problem  of  etablllty  of  laminar  motion  with  respect  to  Infin¬ 
itesimally  small  disturbances*  It  should  be  noted  that 
instability  does  not  neessserlly  lead  to  turbulent  motion*  it 
could  lead  to  another  state  of  laminar  motion*  Another  proULam 
Is*  of  course*  the  turtulent  flow  Itselfo 

The  ordinary  theory  of  stability  of  laminar  flow  has  been 
extended  to  include  the  effects  of  magnetic  forces*  All  of 
these  inyeetlgatione  consider  those  cases  where  the  basic  flow 
is  not  affected  by  the  mantle  forces  but  the  distorbancee^^^* 
Ja  all  these  cases*  it  was  found  that  the  magnetic  field  tends 
to  Increase  the  stability  of  laminar  flow  because  the  dietox^ 
bsneee  tend  to  be  dai^md  out  by  the  eddy  currentso 

It  is  well  known  that  the  stability  of  the  flow  Is  greatly 
affeeted  by  the  profile  of  the  basic  flow*  It  will  be  of  interest 
If  the  effect  of  stability  of  laminar  flow  by  the  magnetie  field 
should  be  investi^ted  when  both  the  effeots  of  induced  naeaetlo 
force  end  of  the  change  of  basic  velocity  profile  occur*  c.g** 


the  problem  of  parallel  flow  with  a  transTeree  magnetle  field 
dlecuaeed  in  ^  lio 

Another  interesting  point  Is  that  the  Squire's  theorenl^ 
fbr  the  relation  of  three  dimensional  and  tw  dimensional  dls* 


turbanees  in  the  stablli^  problem  of  parallel  flow,  does  not 
In  general,  hold  in  magnstohTdrodynamicSe 


d.  Ihrbulsnce  In  Magnetohfdrodynamlcso 

The  interaction  between  the  elsetroeaeaetie  and  the  hgdroo 
dynamic  forees  in  an  Inooi^pressible  conducting  fluid  In  turbulent 
motion  was  first  studied  ^stematically  by  Batcheler^^  and  was 
farther  developed  Ghandraeeld»r^^  and  othersa  the  fundamental 
equations  are  still  equations  (3*1)  to  (3*3)o  Vb  re-rlsw  some  of 
the  rseulte  of  turbiilent  flow  as  foUowai 

(a)  Gbowth^^  of  eleotronagnstle  energy  in  a  turbulent  motione 
Batdielor  considered  the  distui^aixse  of  electromagnetic  energy 
produeed  by  the  turbulent  motion  in  the  absence  of  external  else* 
trloel  or  msgnstlc  fisld*  Ho  assumsd  that  both  ths  tuxbnlsnt 
velooity  and  the  eleetroeegnstic  fields  are  stationary  random 
functions*.  From  equation  (U*13),  it  may  be  shown  that  the  rate 
of  change  of  the  average  amount  of  ma^tio  energy  in  unit  volume 


of  fluid  Is 


A 


(6.1) 


where  the  subeeript  H  refers  to  the  eomponsnt  in  ths  direction 
of  Uie  Unas  of  ma^ietic  force  and  the  repeated  suffix  1  indicates 
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summation  ovtir  three  orthogonal  componentso 

]h  the  absence  of  electronagnetlc  effer'^tsy  for  Isotropic 
tur^ulsncay  It  is  uell  knovn  that  the  mean  square  vorticlty 
satisfies  the  equation 


(g^2) 


Here  the  analogy  of  ma@r»tic  field  and  -vorticlty  may  be  applied 
by  comparing  equations  (6,1)  and  (6.2)«  From  experiments  of 
ordinary  Isotropic  turbilencoy  It  was  found  that  the  rate  of 
change  of  is  approximately  zeroy  i.eo  the  production  term 
^  and  the  dee^  tern  aJ  jvci>/l^  are  ^proxlmately 
in  equilibriumo  If  one  assures  that  the  statistical  distribi»» 
tions  R  are  approximately  the  came  as  those  of  y  the  tso 
contributions  to  d  ('^j/dt  will  be  in  equilibrium  ^  S 

j  |/?^|  will  d^eay  to  aero  and  lfi^^<A)  ,  j  vUl 
inoreaseo  Hence  the  criterion  for  the  energy  of  the  disturbing 
Toasietlc  field  to  increase  is 


>  i  (6,3) 

(b)  The  Inrariant  theory  of  isotropic  turbulence  in 
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magnetohydrodynandcs  • 

Chandrasekhar  extended  the  statistical  theory  of  isotropic 


tuniblenee  of  ordinary  hydroctynamies  developed  by  Tkylory 
von  Karnan  and  others  to  the  case  with  the  presence  of  elaetzo- 
ma^ietie  field.  It  is  convenient  to  vrlte  the  fundamental 
equations  (3.1)  to  (3«h)  in  the  Iblloving  fora 


lAiere  h^  =  ('^')^  Hj.  hj^  has  the  dimensions  of  a  veloeiti7o 
Subscrlnts  i  and  k  refer  to  the  ith  and  Icth  component  of  tbs 
vector  respectively  and  the  summation  convention  is  used*  In 
the  derivation  of  (6ali)  the  solenoldal  properties  of  tbs  vsetors 
7«nd  are  uscido 

3h  the  statistical  theory  of  homogeneous  and  isotropic 
torbulenoe  of  Taylor  and  von  Karman«  we  consider  the  oorrelstlon 
betussn  the  fluctuating  quantities  at  two  points  (pCs^)  and 

In  spacso  Prom  equation  (Sob),  we  see  that  the  follow* 
Ing  correlations  should  be  considered: 

Double  correlations t 

IMple  correlations: 

Hi 

_ _  _  (6.6) 

Where  tbs  primes  and  tbs  laek  of  primes  refer  to  tbs  qaentitlss 
St  tbs  point  end  P  respsetlveljo 


*-‘JL8  = 


For  homo^noous  and  isotropic  turbulence^  mb  nay  apply 
the  thBOiy  of  invariants  to  the  tensors  of  (6»5)  and  (6,6)a 
For  double  eoi*relatlons>  ve  have 


4^,  O' 


^  '  )r  -h  ^  Q)  Sjy, 

U}^  -  h 


(6.7) 


where  z  z  X  f/  *  Q*  H  and  R  axe  scalar 

functions  of  r  and  t  onlyo  t  y  ~  /j  7  ^ 

is  the  usual  alternating  syaboln  j  H 

Sinilarly  for  trlpla  correlations^  ve  have 

-  f  Tt  f/ J  *r^Sn}*i  ts,-^  5 

-  pS  K'  ^  -(fs  tit 


Ui 


^;Jl/  <// 


-  ^c^'-yx-  y  ) 


(6.8 


^ihore  T,  S,U,  V,  P  and  are  scalai's  which  are  functions  of  r 
and  t  in  homogeneous  and  isotropic  turbulence®  Prime  on  these 
scalars  refers  to  partial  differentiation  with  respect  to  r^ 


Prom  the  funda 'cental  equation  {6»U)»  the  equations  gpT- 
emlng  these  scalar  functions  may  be  found*  The  equation  gorsming 


Qg  5g  and  T  is 


c)<p 

c>r 


+•  o 


0(T-^) 


fr)<^ 


(6.9) 


This  la  a  generalization  of  the  ron  Karman  -  Hovar^  equation  to 
nagnetohydrodynanlcs.  Similar  to  ordinary  hydrodynamics ,  equation 
(6.9)  glees  an  invariant 


f  Q(r.J-)r''Jr 

o 


CONiTANT 


(6*10) 


This  Is  known  as  Loltelansky  invariant*  The  existence  of  such  an 
invariant  shows  that  no  transfer  of  energy  from  the  velocity  field 
to  the  magnetic  field  takes  place  among  the  largest  eddies* 

The  equation  governing  H  and  P  is 


(6.U) 


and  the  equation  gsveming  U,  V  and  W  la 

5/? 


FMm  equations  (6.11)  and  (6.12)^  w  nay  obtain  the  invariants  of 
the  loitsian^y  type  for  H  and  R  respeotlvely* 

The  rate  of  dissipation  of  energy  may  be  obtained  txom  (6,9) 
and  (6*10)  taking  the  limit  of  r*~»0.  Vb  have  then 


i^/f7  = 


—  (pc  k)  C) 


(d.i3) 

or 

+-  i'  (6.11i) 

liiero  Q2  and  define  the  curvatures  of  the  curves  of  longitudinal 
correlations  t{,  and  at  r  *  0  respectively.  They  repiasent 

the  snallest  eddies  in  the  turbulent  field.  Equation  (6.11^)  shova 
that  the  rate  of  dissipation  is  due  to  viscous  dissipation  and 
production  of  Joule  heat  by  electrical  conductivity. 

7.  Magnetolgdrodynaiaic  experiments. 

There  are  very  few  expei'iiaental  investigations  of  ngneto* 
hydrodynanic  phenoiaena  because  it  is  necessary  to  have  a  strong 
nagnstic  field  vlthin  a  large  volume  in  order  to  obsenrs  any 
appreciable  magfietol^rdrodynamic  effects.  It  is  rather  difficult 
to  achieve  in  a  laboratory. 

The  first  raagnetohydrodynanic  experiment  was  done  by  Hartmann^* 
7 

and  Lazarus  .  They  investigated  tho  flow  of  msrcuiy  in  a  pipe 
perpendicular  to  a  magnetic  field.  They  found  a  qualitative  agree¬ 
ment  between  theory  {  f  h)  and  experiments  for  laminar  flow  mid 
discovered  that  the  magnetic  field  has  an  Influence  on  the  transition 


between  laminar  and  tiirbulent  flovo 

Lehnert^  recently  reaxaminec!  Hartmann's  data  and  fouixi 
that  the  critical  velocity  for  transit  ion  was  proportional  to 
the  magnetic  field  strengths  Lohnert  also  measured  the  torque 
transmitted  by  mei*cury  contained  between  two  non-conducting 
cylinders  rotating  in  a  magnutlc  field  parallel  to  the  axiso 
For  laminar  flow,  the  effect  due  to  magnetic  field  is  smalle 
For  non-laminar  flow  with  the  outer  cylinder  at  rest  and  Inxer 
cylinder  rotating,  the  to]*quo  was  found  to  decrease  when  the 
mastic  field  was  applying,  lehnert  planned  to  use  molten- 
sodium  instead  of  mercury  in  his  new  experiment  because  of 
higher  conductivity  and  lower  density  than  mercury.  For 
instance  the  value  of  for  molten  sodium  is  about  35  time* 

that  of  mercury  in  the  same  ptysical  conditions.  However, 
there  are  many  experimental  difficulties  to  use  molten  sodium. 

Lehnsrt^  also  showed  a  very  interesting  simple  s^qwrinsnt 
to  demonstrate  the  mag^etohydrodynamlc  phenonsnon.  Consider  a 
cylindrical  glass  vessel  containing  mercury.  Without  magnetic 
field,  if  one  moves  the  amreury  with  a  peg,  the  surface  is 
agitated  and  complicated  wave  pattern  occurs,  particularly  vhen 
one  hite  the  vejE<sel,  surface  waves  occur  Just  like  water  surlhce 
waves. 

If  the  vessel  is  placed  in  a  strong  magnetic  field  perpen¬ 
dicular  to  the  surface  of  the  liquid,  the  surface  waves  disajqjeer. 
The  liquid  no  longer  behaves  like  water  but  like  thick  syrup. 

If  a  peg  is  being  moved  in  the  liquid,  the  surface  shows  some 


largo  lAirls  all  with  an  axis  of  ixjtatlon  parallel  with  the 
magnetic  field  llne^. 


Part  IL  MAG’ClOGASDYrZAr'flCS 


At  very  hl^  temperature,  the  gas  nay  be  ionized,  and  the 
Interaction  of  electromagnetic  forces  and  the  gasdynaraie  force 
may  not  be  negligible*  Purthemiore,  under  such  conditions,  the 
compreeelblllty  effects  should  be  considered.  In  this  ease  our 
fund^nental  equations  are  (2.1)  and  (2.9).  This  system  of 
equations  has  not  been  Investigated  yet..  Only  the  siaqple  one- 
dlannslonal  waves  of  Infinitesimal  amplitude,  Alfven's  wave 
(  ^9)  and  shock  waves  (4 10)  have  been  discussed  under  further 
slnqilifled  assumptions,  kh  shall  briefly  review  these  results 
later.  First  we  should  like  to  bring  out  the  important  para¬ 
meters  in  ma^ietogasdynanics.  We  shall  restrict  ourselves  to 
ths  adiabatic  ease  s  o)  with  negligible  radiation  lost 

(IriS  0).  Fuidhermore,  we  shall  consider  the  case  where  the 
velocity  of  the  flow  is  much  smaller  than  the  velocity  of 
Uf^t  C  80  that  the  relativistic  effects  may  be  ne^ected  and 
that  the  energy  in  the  eleetrle  field  is  negllgiblB  in  congArlaon 
with  that  in  the  magnetic  field. 

Ohder  these  nagnetogasdynamic  approximations,  the  Ihndamental 


equations  (2.1)  to  (2.9)  may  be  reduced  to  the  following  non- 
dimensional  formx 


(8.2) 


^r 


(8.3) 


^/v)V^  =  r  T 


(8.1i) 


n^'Rrf[^*^-(^"‘f')J^(r-0 


h 


J»» 


«h«  r*.  I  I  T»  =  ,  i;"*.  (>. 


J  '^>/ “  i/'«  ^  ^^ty' 
hf*-Jl  ys.-^  ,  -5 _ ^ 


(8.6) 


Mach  naiher  s  m  *  — I.  »  *  Rrandtl  nuid)er.  and  the 

a^/ITi  ^  f\o 

other  non<^ln8n8lonal  quantltloB  are  saaa  as  thoee  defined  in 
(3.U). 

From  these  equations  «e  see  tiiat  as  far  as  the  magoetlo 


forces  are  concerned^  the  Important  parameters  are  still  ^  and 
^  The  other  new  parameters  are  those  well  known  parameters  la 
gas  dynamics,  i.e.  P,.,  M  and  ^ . 


The  general  solution  of  magneto gasdynaialcs  is  vary  diffi* 
cult.  In  order  to  bring  out  some  essential  fe  itures  of  magneto* 
gasd^naralcsf  we  may  investigate  some  very  single  cases*  One  of 


these  simple  cases  is  one  dimensional  wave  motion  of  small 


aiqilltude*  Alfven  was  the  first  one  who  found  such  wave  in 
magnetohydrodynamlcs'''”*  ^  there  is  a  homogeneous  magnetic 
field  in  an  Incompressible  and  Inviscld  fluid  of  density  ^ 
and  infinite  conductivity^/" ,  the  disturbance  in  this  liquid 


will  propagate  as  a  wave  in  the  direction  of  with  a  speed  of 


Ihls  wave  is  known  as  Alfven's  waveo 

Alfven *8  analysis  has  been  extended  to  the  case  of  cooqjres- 

sible  fluid  by  several  authors  *  The  following  analysis  was 
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dua  to  van  de  Hulst  . 

Strictly  speaking,  to  investigate  the  one  dimensional  wave 
motion  in  magnetogasdynamies,  we  should  use  the  system  of  equation 
(8*1)  to  (8*^)*  But  Van  de  Hulst  made  the  following  simplified 
assumptions  X 

(1)  The  energy  equation  is  replaced  by  a  sinpis  adiabatic 


relation 


'= -Cd.' c/ 

where  is  the  perturbed  pressure  from  that  of  the  gas  at  rest 
^  ,  IV  is  the  material  displacement  from  the  position  at  rset 


and  the  fluid  velocity  q  is 


2 

C  s  ”  compressibility  factor,  is  the  velocity  of  sound 

of  the  gas  at  teste 

(11)  The  ccefflcient  of  viscosity  is  assumed  to  be  eonstanto 
(iii)  All  the  quantities  are  assumed  to  be  function  of  one 
space  coordinate  2  and  time  t  only.  The  perturbed  quantities  are 
assumed  to  be  small  so  that  the  equations  for  the  perturbed  quan» 
titles  may  be  linearized. 

(It)  a  li^geneous  magnetic  l^sld  H  is  ^plied.  Itlis  wish  to 
investigate  wave  travels  in  the  Z-dlrectlon.  Tne  x««xls  is  <dM>sen 
perpendicular  to  H  so  that  the  components  of  H  are  (0, 

If  we  substitute  the  assumptions  into  equations  (2.1),  (2,2), 
(2*3 (2.6),  (9.2)  and  (9«3},  simplify  them  under  magnetogas* 
dynamic  c^proxlmations  and  linearized,  we  have  two  independent 
sets  of  linear  equations  for  the  perturbed  quantities;  one  set 
contains  the  perturbed  quantities,  hy»  Ex>  ^y* 

9  and  the  other  set  contains  h^,  Ey,  «^,  q^,  and 
Vfe  are  looking  for  periodic  solutions,  in  uhich  all  perturbed 
quantities  are  ifuroportlonal  to 

where  ^  has  a  given  x«al  value.  Only  for  certain  eigenvalues 
of  k,  such  solution  exists.  These  eigenvalues  of  k  give  the 
different  modes  of  wave  propagation  and  they  are  detemined  from 
the  deteimlnantal  equations  of  the  two  sets  of  linear  equations 


for  the  perturbed  quantltieso  Ihese  two  detemlnantal  equations 
after  simplifications  are  aa  follows 


iriiare 


(9.6) 


(9.7) 


Since  -f  and  f  ^are  linear  functions  of  x,  equation  (9.$)  has  three 
roots  of  X  and  (9.6)  has  two.  These  roots  are  the  eigemraluss  we 
are  looking  for.  Ws  shall  not  review  all  the  special  cases  for  z 
here  but  a  few  Interesting  ones* 

(a)  Ifo  external  aagnetie  field. 

The  roots  of  (9*5)  are 

*  *  ^  ^  '■»  *  damped  eleetronagnetie  wave  (9*8a) 

z  s  z^  s  ^  - ^  •  viscous  wave  (9*8b} 

z  s  ^  dangled  sound  ware  (9*6o) 

There  is  x»  coupling  between  these  Uiree  fundamental  modes*  The 
roots  of  (9*b)  ate  those  of  (9*8a)  and  (9*8b}*  The  differenoe  of 
nodes  with  the  same  values  of  z  is  a  difference  of  the  plane  of 
polarisation* 
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(b)  Propagation  In  tbs  direction  of  the  external  field 

g  *  Oe  The  soiuid  waves  are  not  affected.  The  other  modes  are 
coupled  and  follow  from  the  equation 

(c)  Propagation  perpendicular  to  the  external  field. 

1  «  0  but  g/x^^O«  The  viscosity  wave  of  (9.5)  and  the  electro- 
megnetie  wave  of  (9.6)  are  not  affected.  The  viscosity  wave  of 
(9.6)  la  slightly  affected.  The  two  other  isqjortant  modes  of 
(9.5)  are  coupled  and  satisfy  the  equation 


)  "*  ^  (9*10) 

(d)  Ihidamped  waves.  If  the  conductivity  is  aero,  we  have  the 
ordinary  undamped  light  and  sound  waves. 

If  the  damping  terns  are  omitted,  which  correspond  to  infi¬ 
nite  conductivity  ( d"  •=>  «>o  )  end  zero  viscosity  ^  :  0),  the 
equations  (9.^)  and  (9.6)  reduce  respectively  to  the  following 
simple  form 


2±. 

x-x^ 


(9.U) 


and 

X  ~  ^  -^0  X<.  J  (9.12) 

The  two  solutions  of  (9.11)  depend  on  the  relative  magnitude  of 
sound  T^oeity  of  the  nsdium  a^  and  the  veloolty  of  tha  magnato- 
hydrodTsamie  wave  V  There  is  a  slow  mode  with  velocity  smaller 


than  sound  veloolty  and  a  fast  mode  with  velocity  larger  than 
sound  velocityo  If  the  magnetic  field  increases^  the  slow  mod* 
changes  from  the  nolution  x  a  Xj^  to  the  solution  x  *  (g  1) 

Ihis  mode  represents  a  regular  Alfven's  wave  with  velocity  Vg  as 
long  as  Vjj  is  much  smaller  than  a^»  It  changes  to  a  kind  of 
regarded  sound  wave  fbr  larger  values  of  The  fast  node  starts 
o  t  as  an  ordinary  sound  wave  but  changes  to  a  modifiad  Alfrsn's 
wave  as  soon  as  Is  larger  than  a^» 

(e)  Slightly  damped  vave«  If  both  "a"  and  "b"  are  snallf 
gone  first  order  effects  due  to  damping  nay  be  found  from  (9r5) 
and  (9«6)o  Ve  will  not  discuss  them  hereo 

10«  Shock  waves. 

In  the  last  section^  we  consider  the  propagation 

of  waves  of  infinitesimal  amplitude  in  a  conqpiressibXB  fluid.  The 

sttidy  of  waves  of  finite  amplitude  in  nagnetogasdynamlcs  has  not 

been  carefully  carried  out.  Only  the  analysis  of  analogues  of 

the  Rankine^Hugonlat  equations  for  shock  waves  in  an  infinitely 
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conducting  fluid  has  been  mads  ^  de  Hoffman  and  Teller  «  Ihey 
discussed  both  ralatlvistic  and  non^relativlatic  shock  vanres* 

But  ve  shall  restrict  ourselves  to  ths  non-relatlvlstic  cases 
which  have  been  systematically  interpreted  by  Helfer^^, 

One  of  the  main  assumptions  in  de  Hoffman-Teller  analysis 
is  that  the  conductivity  of  the  fluid  is  infinite.  Even  thou^ 
the  conductivity  of  the  fluid  is  actually  finite  because  of  the 


large  dimension  of  Interstellar  clouds  In  which  such  an  analyals 


is  api^led  and  the  large  conductivities  of  stellar  materials^ 
such  an  assun^tion  is  a  reasonable  onoe  Under  this  assumption 
the  lines  of  magnetic  force  move  with  the  fluid  U(e^« 
Another  assumption  is  that  the  fluid  le  aesisiied  to  be  an  ideal 
gae  of  a  constant  value  of  ratio  of  specific  heats  'iT  • 

Fbr  ordinary  shock  relations »  if  we  know  the  conditions  in 
front  of  the  shocks  the  conditions  behind  the  shock  depend  on  one 
paramster,  the  strength  of  the  shock,  e.g.  the  pressure  ratio 
across  the  shock.  In  the  case  of  magneto gasdynamics,  the  shock 
relatione  depend  on  three  parameters: 

(l)  Strength  of  the  shock.  We  may  use  the  pressure  ratio 

Y  *  Pz/pl  Aa  a  measure  of  the  strength  of  shock. 

(ii)  Migpetlc  field  strength.  Vo  may  use  the  ratio  of 
magnetic  energy  per  unit  mass  to  internal  energy  Q(^  -1)  as  e 
msssure  of  the  importance  of  the  magnetic  field  strength,  ihers 


Q  =  r~  (ID 

(ill)  The  direction  of  the  magnetic  field.  We  mi^  use  the 


(10.1) 


angle  between  magnetic  Unas  of  foroa  in  front  of  the  ahoek  aatf 
the  direction  of  propagation  of  shock  &  or 

83^  *  tan  d  1  (10,2) 

as  a  maasura  of  the  direction  of  magnetic  field.  Ihe  sUbaexipt  1 
refera  to  vdLuee  in  float  of  shock  and  svbaerlpt  2  rafers  to 
valOM  behind  the  shodr. 


^30- 


According  to  the  paraneter  we  have  the  following  three 
cases t 

(a)  Parallel  shock  Sj  S  0 

]h  this  case*  tho  magiMtic  field  has  no  influence  upon  the 
gas  4jmamie  phenoni9na«  Classical  Rankine-Hugoniot  relations  for 
shock  hold  in  this  casoo 

(b)  Perpendicular  shock  e^ 

In  the  coordinate  system  of  stationary  shock  with  velocity 
Qj  In  front  of  the  shock  and  perpendicular  to  the  shock  front,  we 
have  the  following  shock  relations  in  magneto gasdynandca* 

^>0,  (10,3) 


(10.1*) 


f<  % t^Hr.  r  ^ 


K 


'  ^  -  X'  + 

^  AT/ 


(10.6) 


(W.7) 

Equatione  (10.3)  and  (10.5)  have  exactly  the  same  forme  as  la  tte 
case  of  ordinary  shock  if  one  takes  into  account  the  mgnetie 
pressure  and  energy  by  (10.6)  and  the  coupling  of  magmtic  and 


velocity  field  by  (10.7)  whl(ih  is  a  consequence  of  infinite 
conductivity.  The  magnetic  field  is  compressed  by  the  shock 
exactly  to  the  same  extent  as  the  fluid 


Ihe  equation  for  in  terns  of  T  and  Q  is  obtained  frma 
(10,3)  to  (10.7)  as  followei 

When  <3  s  0,  it  reduces  to  the  Rankine-Rugoniot  relation.  The 
presence  of  the  magnetic  field  is  seen  to  cause  a  decrease  in 
compression  ^  .  The  Mach  numbers  q2/*j^  and  q2/a2  increased 
due  to  the  presence  of  the  magnetic  field  for  the  sans  reason, 

(o)  Oblique  shocks 

The  oblique  shock  relations  of  this  case  has  been  worked 

out  by  Heifer.  NUnerieal  results  for  Q  vs  T-1  at  various  values 

of  s^  are  given  by  Heifer  in  reference  2Uo  It  is  interesting  to 

note  that  for  very  weak  shock,  t'^1,  the  shock  relations  tend  to 

be  the  value  of  wave  of  Infinitesimal  amplitude  discussed  in  $  9* 

Another  Interesting  result  is  that  for  weak  fields,  the  magnstle 

field  is  always  amplified  by  the  passage  of  shock  fronts.  This 

23 

Is  important  to  explain  some  of  the  astroi^iyslcal  phenomena 
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